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ABSTRACT 

K band observations of the galaxy populations of three high redshift (z = 0.8-1.0), 
X-ray selected, massive clusters are presented. The observations reach a depth of 
K ~ 21.5, corresponding to K* + 3.5 mag. The evolution of the galaxy properties 
are discussed in terms of their K band luminosity functions and the K band Hubble 
diagram of brightest cluster galaxies. 

The bulk of the galaxy luminosities, as characterised by the parameter K* from 
the ISchechterl l)1976j) function, are found to be consistent with passive evolution with 
a redshift of formation of Zf w 1.5-2. This is consistent with observations of other 
high redshift clusters, but may be in disagreement with galaxies in the field at similar 
redshifts. A good match to the shape of the Coma cluster luminosity function is found 
by simply dimming the high redshift luminosity function by an amount consistent 
with passive evolution. The evolution of the cumulative fraction of K band light as a 
function of luminosity shows no evidence of merger activity in the brighter galaxies. 

The evolution of the brightest cluster galaxies (BCGs) is tested by their K band 
Hubble diagram and by the fraction of K band cluster light in the BCGs. The evolution 
observed is consistent with recent previous observations although the scatter in the 
Hubble diagram allows for a range of evolutionary histories. The fraction of cluster 
light contained in the BCGs is not smaller than that in Coma, suggesting that they are 
already very massive with no need to hypothesise significant mergers in their futures. 
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1 INTRODUCTION. 

The evolutionary history of galaxies in clusters remains a 
subject of debate. The two most common explanations of 
massive early-type galaxy formation and evol ution are those 
of monolithic co llapse (e.g. l Eggen et alll962l) and hierarchi- 
cal merging (e.g. lCole et al.ll994l) . In the monolithic collapse 
scenario all the galaxies (and the stars therein) are formed 
in a single burst and subsequently evolve passively (along 
the main sequence) with no further star formation. Such 
a model will result in a very homogeneous population of 
galaxies since their ages and metallicities will all be close to 
identical, with a small degree of scatter reflecting the varia- 
tions in initial mass function at formation. There is a signifi- 
cant body of observational evidence that luminous early type 
cluster galaxies are indeed remarkably homogeneous. For ex- 
ample, the observed tightness of the colour-magnitude rela- 
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tion (e.g. lVisvanathan fc Sandagdll977l . iBower et allll993l 

is naturally explained by such a homogeneous population. 

In the merging model galaxies form by a series of merg- 
ers within a hierarchic al model of structure formation (e.g. 
Kauffmann et al Jl995tl . The hierarchical scenario presents a 
radically different evolution than in the monolithic collapse 
picture. Bursts of star formation, related to mergers, may 
occur and a more gradual increase in the number of mas- 
sive galaxies in a cluster would then occur as small galaxies 
merge to form larger ones. 

Although the ages of the majority of the stars in galax- 
ies in both scenarios are similar (in order to explain the 
tightness of the colour-magnitude relation and the evolu- 
tion of the fundamental plane), the mass as a function of 
look-back time will be quite different, i.e. the epoch of the 
assembly of massive galaxies is not the same as the epoch of 
star formation in the merging case. The hierarchical model 
will exhibit negative evolution of mass as a function of red- 
shift (i.e. at higher redshifts the massive galaxies will be 
less massive on average), whereas in the monolithic collapse 
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model galaxies will have a constant mass with redshift. The 
negative evolution expected from hierarchical models will 
be most apparent in the most massive galaxies, which are 
predicted to have assembled more recently. The redshift of 
assembly of massive galaxies is, however, poorly known and 
is probably dependent on environment. 

If one is interested in the stellar mass of the galax- 
ies, as opposed to their star formation rate, then the K 
band is a good choice in which to make observations since 
the light at such wavelengths originates mainly from the 
longer lived stars of the ma in sequence (see e.g. figure 1 of 
Kau ffmann fc Charlodll99^) . and has the added advantages 
that the k-correction differences between galaxies of different 
spectral types are small in K and the Galactic extinction by 
dust is small. Thus computing K band luminosity functions 
for clusters of galaxies should give a useful measure of the 
mas s distribution of galaxies wit hin clusters (see section 2] 
and iTrentham fc Mobashe'rlll998l) . 

Brightest cluster galaxies (BCGs) are known to have a 
very limited variation in absolute magnitude which histor- 
ically led to their nomination as a candidate for a cosmo- 
logical standard candle wi th which to directly measure th e 
curvature of the Universe lSandagelll972al . lSandagdll9 72b.') 
Controversy over their suitability as standard candles is be- 
lieved to be due to environment, with BCGs in high Lx clus- 
ters e xhibiting much less scatter than those in low Lx sys- 
tems ferough et al.ll200l iBurke et ai1l2000ft . However the 
evolution of BCG properties with redshift is now the pri- 
mary interest of research on BCGs, since there is much evi- 
dence that they are a special case in the evolution of galaxies 
within clusters. For instance, it is known that BCGs often do 
not follo w the same luminosity function as other galaxies in 
clust ers iSchechtedll976l . |PressleJll978l iBhavsar fc Barrowl 
1985), most probably due to the fact that they have a pe- 
culiar formation history. Therefore the evolution of BCGs 
can provide a different and complementary study of galaxy 
formation theories compared with the general cluster pop- 
ulation. The K band Hubble diagram of BCGs provides an 
efficient measure of the evolution of BCGs, and we present 
our results of this in section [5] 

Throughout this paper we have used a cosmology of 
Ho = 70 km s" 1 Mpc -1 in a flat universe with S1m = 0.3 
and f?A = 0.7, except where stated. 

2 DATA. 

2.1 The Clusters. 

We present a study of three of the most mass ive (~ 1O 15 M0 
iMaughan et al J 12001 iMaughan et al J 120031) . high redshift 
clusters known. They are thus ideal probes of galaxy evolu- 
tion within clusters. Such massive clusters at high redshift 
are very rare and we have an unusual opportunity to study 
the galaxy populations of rich, distant clusters and com- 
pare results with local massive clusters such as Coma. The 
high redshift of the clusters should make any evolution in 
the galaxy populations easier to observe. Two of the clus- 
ters (C1J1226 and C1J1415) appear relaxed based on their 
X-ray morphologies and one of them (C1J0152) is probably 
in a state of merging. Thus we also have a small selection of 
different environments. All three clusters are X-ray selected 



and thus should be relatively free from projection effects. 
The galaxy properties should also be largely independent of 
the selection process. They were all disc overed in the Wide 
Angle ROSAT Pointed Survey (WAR PS jScharf et alll997t 
Jones et al. 1998, Perlman et al. 2002). Details are given be- 
low on the X-ray properties of each cluster. 

2.1.1 CI J 0152 

The cluster C1J01 52 7-1357 jEbeling et al] |200A 
iDella Ceca et al] l2000t IMaughan et al] 120031) is at a 

redshift of z — 0.833 and has a bolometric X-ray luminosity 
of 1.6 ± 0.2 x 10 45 ergs s _1 . It is composed of two major 
subclumps which ar e probably gravitationally bound 
iMaughan et al] l2003Tl . The projected separation of the 
subclumps is ~1.5 arcmin (or 720 kpc), suggesting that they 
may be in the early stage of a very massive merger event. 
Individually the subclumps have X-ray temperatures and 
luminosities of L x = 1.0 ± 0.2 x 10 45 ergs s _1 , T x = 5.5±q J 
keV and L x = 5.8tJ;£ x 10 44 ergs s~\ T x = 5.2± ; 4 keV 
respectively. Thus, each clump considered individually is 
still a massive cluster with M tot ai ~ 6 x 10 14 Mq and 
Aftotai ~ 5 x 10 1 4 Mp) for the northern and southern clumps 
respectively (see IMaughan et al]l2003ll . 

2.1.2 CI J 1226 

C1J1226.9+3332 jEbeling et al.ll200ll IMaughan et"aTll200l 
ICagnoni etafl 1200 J) is at a redshift z = 0.888. It has a 
bolometric X-ray luminosity of Lx = 5.3^ 2 x 10 45 ergs 
s" 1 , an X-ray temperature of Tx = 11.5t.n'a keV and 
Mtotai = 1.4loia x 10 15 M Q iMaughan et al]l2003l) . making 
it one of the most massive clusters known at high redshifts. 
The cluster appears remarkably relaxed in X-rays indicating 
little dynamic activity. The K band image with adaptively 
smoothed XMM X-ray contours (from Maughan et al 2003b) 
is shown in figure 

2.1.3 CIJ1415 

C1J1415. 1+3612 is at t he highest redshift of the three clus- 
ters with z = 1.03 JPerlman et al] l2002h . Although at 
present we have only ROSAT X-ray data for the cluster 
so sub-structure may not be resolved, it exhibits smooth X- 
ray contours suggesting it may be dynamically relaxed. It 
has a bolometric X-ray luminosity of L x = 2 x 10 45 ergs 
s _1 , and is thus likely to be a very massive cluster. 

2.2 Observations. 

All three clusters were observed at the 3.8m United King- 
dom Infrared Telescope (UKIRT) on Mauna Kea, Hawaii 
using the UFTI camera, a 1024 x 1024 pixel HgCdTe ar- 
ray, with a pixel size of 0.091 arcseconds. Observations were 
made with the K98 and J98 filters (50% cut-offs 1.17-1.33/xm 
and 2. 03-2. 37pm respectively) with an exposure time of 66s 
per frame (or 60s per frame for a minority of observations 
made in service time) and the seeing was typically w 0.5 
arcseconds in K. The total exposure times in K are listed in 
table^along with the average seeing for each cluster. A spa- 
tial dithering of 15 arcsec was employed. Since the field of 
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Figure 1. UKIRT UFTI K band image of the z = 0.89 cluster 
C1J1226 with adaptively smoothed XMM X-ray contours from 
Maughan ct al. (2003b). The image is ~ 2.9 arcmin square. 



view of 90x90 arcsec covered only a fraction of each cluster, 
mosaics of different si zes for each cluster wer e used. On each 
night standard stars iHawarden et al.ll200ll) were observed 
in order to calibrate the zero-point, air mass coefficient, and 
colour coefficient. Some nights had patchy cloud cover, and 
thus yielded inaccurate photometric calibration. Such nights 
were deemed to be non-photometric, and a calibration was 
obtained via observations of the same fields on photometric 
nights. Non-photometric data were then weighted according 
to their shorter effective exposure times. 

Observations were also made of one or more offset fields 
in the vicinity of each cluster on the sky, chosen as a reason- 
ably 'blank' part of a digitised sky survey image (since at 
these redshifts, the clusters are also blank on digitised sky 
survey images). The purpose of these fields was to provide 
a means of estimating the contamination from foreground 
and background galaxies when constructing luminosity func- 
tions. The offset fields were typically 4 arcmin (1.9 Mpc) 
from the clusters, beyond the virial radii. 

Observations of C1J0152 were made in service time, on 
the nights of 15th and 20th October 2000, 24th Decem- 
ber 2000, 28th December 2000 and 16th January 2001. The 
northern and southern clumps will be referred to as field A 
and field B. The nights of 16th January, 24th December and 
28th December were found to be non-photometric and were 
treated accordingly. The areas observed were 1'50" x 1'39" 
and 1'51" x 1'50" for fields A and B respectively, correspond- 
ing to a physical size of 836 x 752kpc and 843 x 836 kpc. 

C1J1226.9+332 was observed on the 18th, 19th and 20th 
April 2001. The night of the 18th had patchy cloud cover 
and was deemed non-photometric. The observations covered 
an area of 2' 47" x 3'l", corresponding to a physical size of 
1.3 x 1.4 Mpc. 

C11415 was observed on the 10th, 11th, 12th April 2002. 
All nights were photometric with an average seeing in K of 



0.52". The observations covered an area of 2'2" x 2'44", 
corresponding to an area of 1.0 x 1.3 Mpc. For this cluster, 
with the highest redshift, two offset fields were observed in 
order to reduce the uncertainties in the field subtraction. 
For the other clusters, one offset field was observed. 

2.3 Data Reduction. 

Each observation of a particular field consisted of a num- 
ber of dithered frames of the object plus a dark frame con- 
taining only the instrumental background. The data were 
dark-subtracted, flat-fielded and combined using standard 
iraf procedures. A bad-pixel mask was applied to exclude 
any defective pixels. Flat-field frames were produced from 
the object frames, using the median value of each pixel after 
applying a 4a clipping and scaling each frame by its modal 
value to account for sky brightness variations. 

Each frame had its background sky value subtracted, 
using a 3<r-clipped mean as an estimate of the sky value. 
The relative offsets of each frame were then measured using 
bright stars in the field of view. Finally the frames were 
shifted and combined using the mean value of each pixel, 
again after applying a 4<r clipping around the median value. 



3 PHOTOMETRY. 
3.1 Standard stars. 

The photometry of the galaxies was calibrated using stan- 
dard stars observed on the same nights as the cluster fields. 
The true magnitude of any celestial body is taken to be 

m truo = ZP- 2.51og(C7?) + Asecz + B{J ~ K) trua (1) 

where ZP is the zero-point magnitude, CR is the count 
rate, A is the coefficient of extinction per unit airmass, z is 
the zenith angle (secz therefore being the airmass) and B 
is the colour coefficient necessary due to differences in the 
response of the combination of camera and filters used and 
the standard star system. 

Values of ZP, A and B were determined from obser- 
vations of standard stars taken on each night. Standard 
stars were selected fr om the UKIRT faint standards list 
(Hawardcn ct al. 2001). Typically each night 6 stars of vary- 
ing colours were observed between 1 and 3 times each at 
differing airmasses. 

Whilst the value of A and ZP may vary from night to 
night, the value of the colour coefficient B should be almost 
constant. Thus a single value of B was measured from all 
the photometric nights combined. The value we obtained 
is consistent with zero , as expected from the filter design 
jTokunaea et al]|2002l) . An independent value of ZP was 
measured for each night, and for most nights, an indepen- 
dent value of A was measured. Where this was not possible 
a value of A from other nights was used and checked for 
consistency. All were found to be consistent. In any case the 
low extinction in the K band and the low airmasses at which 
most observations were performed combine to make this a 
relatively small correction (<0.2 mag). 

The parameters are listed in tableQ Note that the final 
photometric error, as estimated from the scatter in the stan- 
dard star measurements, is less than 0.1 mag. Corrections 
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Table 1. Calibration parameters for K band. 

C1J0152A C1J0152B C1J1226 C1J1415 



A -0.150 ±0.110 -0.150 ±0.110 -0.178 ±0.052 -0.088 
B -0.0043 ± 0.094 

ZP 22.248±0.018 22.248±0.018 22.367 ± 0.008 22.422±0.044 

Mean secz 1.445 1.403 1.210 1.160 

Exposure/ s 3882 3942 23946 17820 

Seeing/ arcsec 0.53 0.53 0.57 0.52 



due to Galactic extinction are negligible («0.007 mag) and 
were not applied. 



3.2 Galaxies. 

The SExtractor software of lBertin fc Arnouts! (I1996T) was 
used to search for objects in each field. To detect objects 
a threshold value per pixel must be chosen along with a 
minimum number of connected pixels. Because the final im- 
ages are constructed from a jittered pattern of images, the 
depth of observation varied across the final image, being 
at its deepest in the centres and shallowest at the edges. 
Therefore detection of objects was done in two regions for 
each image, using different detection parameters, in order to 
reach as deep as possible in the centre of the image whilst 
avoiding spurious detections at the edges of the image. Typ- 
ically the minimum significance of the detections were 3.5a 
or 4<t in the centres of the images and 7a at the edges, where 
a is the background RMS determined from counts over the 
whole image (and is thus an overestimate of a in the central 
region and an underestimate in the outer region, although 
this is not the case for the central region of C1J1415, when 
a was determined from the central region alone). Reliabil- 
ity of the object detection, and in particular the handling 
of overlapping objects, was checked by eye. Faint pixels sur- 
rounding deblended objects are assigned to one of the sub- 
objects with a probability based on the expected contribu- 
tion at that pixel from each of the deblended objects (see 
iBertin fc ArnoutsllT996l) . 

Counts were me asured using an adaptive aperture based 
on Kron's algorithm jKronl 1980T) , and also in a circular aper- 
ture with radius chosen to maximise the signal to noise. 
For each object a value of stellarity was also measured (see 
IBertin fc Amout Jll993) . 

To determine the magnitudes of the objects equation 
was used with the following complication. It is unknown to 
start with what the true colours of the objects are, therefore 
the magnitudes cannot be determined. To circumvent this 
problem an approximate magnitude was measured in the 
J and K bands neglecting the colour term in equation 
This allowed an approximate colour to be determined. The 
approximate colour is then multiplied by a correction factor 
previously determined using the same technique to measure 
the approximate colours of standard stars and their true 
colours. The average value of , l J "^^ truc was 0.982. 

^ (,J — JV Japprox 

Note that for the determination of colours, magnitudes 
derived from fixed, circular apertures were used, whereas the 
adaptive aperture magnitudes were used to derive pseudo- 
total magnitudes. The use of adaptive aperture magnitudes 
using S Extractor's Kron radius should avoid the prob- 



lem described bv lAndreonl i2002t) of underestimating fluxes 
for galaxies with low central surface brightness. The over- 
all reliability of the photometry was checked by comparing 
the field galaxy counts derived from the offset fields with 
deeper published results (see fig- El • There is generally good 
agreement down to our limiting magnitudes (see below). 

Star-galaxy discrimination was determined using SEx- 
tractor'S stellarity parameter. A cut-off of 0.8 was selected 
to delimit the two classes of objects with those objects with 
values greater than 0.8 being excluded as stars. This value 
was confirmed by examination of the radial profiles of de- 
tected objects. It was found that objects with a stellarity 
greater than 0.8 had an almost constant Gaussian FWHM 
close to the value of the seeing, whereas objects with stellar- 
ity less than 0.8 had more extended profiles. The results are 
insensitive to the precise value of this cut-off as there were 
very few stars in each field. 



4 LUMINOSITY FUNCTIONS. 

K band luminosity functions (KLFs) were calculated sepa- 
rately for each cluster of galaxies. This was done in the fol- 
lowing way. Firstly, some of the brightest galaxies have been 
spectroscopically confirmed as cluster members, so where 
there are no other bright galaxies in the field, no subtraction 
is necessary. At fainter magnitudes, the number of galaxies, 
N c \, in a bin of width AM was counted, with AM = 1 for 
C1J0152 and C1J1415 and AM = 0.75 for C1J1226. The con- 
tamination from foreground and background galaxies was 
estimated by counting the number of galaxies, A^back in the 
offset field(s) associated with the cluster. Thus the number 
of galaxies in each bin is given by, 

iVbin = N cl - AW x (-^-) (2) 

V Aback / 

where A c \ is the area of the cluster field and Ab ac k is the 
area of the offset field. Note that because two different de- 
tection thresholds have been used for each field the faintest 
bins are from a smaller area than the brighter bins (since 
only a small area was deep enough to detect objects of this 
magnitude). Therefore the detections from the smaller area 
have been scaled by an appropriate factor to bring them into 
line with the detections from the total area. 

The corresponding error on the number of galaxies in 
each bin is found by summing in quadrature the error on 
N c i and the error on Ab ac k, 

^Nbin = V erg, cl + CTNback ( 3 ) 
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where cr Nc i is the Poissonian error (lGehrelsl ll98(f). The 
background error, CNback, has a Poissonian term and a term 
taking into account the clustering properties of galaxies, viz., 



Number Counts of Field Galaxies 



ONback 



back X 



2nNA^6't 
2-5 



(4) 



where N is the number density of galaxies in the bin 
and 6 C is the angular radius such that Q = tyO^ where fi is 
the solid angle of the background field. The parameters 
and 8 describe the angular correlation function of galaxies 
such that, 

u>(6)=A u 8 s . (5) 
l|Peebleslll98ol pg.175). 

The value of 8 — —0.8 was used and was found using 
the formula 



log 10 A„ = 7.677 - 0.3297n in 



(6) 



derived by fitt ing to t he data in tab le 1 of 

iBrainerd et alJ (Il995l) . Following lDiorgovski et alJ (Il995l) we 

used nim - K\i m ~ 3. 

Our limiting K magnitude was determined by plot- 
ting the number density of galaxies from the offset fields 
as a function of magnitude and determining the drop-off 
point when compared to large area field surveys. Figure [5] 
shows the number density of galaxies in our four offset 
fields. C1J1226 and C1J1415 were found to be complete to 
isTiim = 21.5 and C1J0152 was found to be complete to 
Klim = 20.0. 

No corrections have been made for gravitational lensing 
of the background galaxies by the cluster dark matter. This 
effect magnifies the fluxes of background galaxies. Trentham 
& Mobasher (1998) show that the net increase in the number 
of galaxies to be subtracted is small (<5%) at the redshifts 
and magnitudes considered here (and for an equally massive 
cluster), because a large fraction of the field galaxies along 
the line of sight are in the foreground. This correction is 
significantly less than the statistical errors. 

Luminosity functions for the three clusters are shown 
in figure |5J 0]and |5] along with the best fitting Schechter 
functions as described below. 



4.1 Fitting Schechter functions. 

The parametric ISchechteil dl97d) luminosity function is of 
the form, 



gdL = /(£re-<*>d(L/L-) 



(7) 



where L* is the characteristic luminosity at which the 
function turns over to the faint end which has slope a. The 

normalisation of the function is given by (j>* . The corre- 
sponding function in terms of K band magnitudes is 



d ^no -^-^ +1 e - 10 °- 4(K *- K: 



AK 



(8) 



where K* is the characteristic K band magnitude. Note 

that <j>* in equation|H|is not identical to the <f>* in equation[7| 
but in fact differs by a factor — 0.41nl0. 

In general the Schechter function is found to be a good 
fit to composite cluster galaxy LFs, as in the 2dF Galaxy 



i l i ® 



Figure 2. Number density of field counts as a function of 
magnitude. Open circles are from the combined offset fields of 
C1J1226 and C1J1415, op en squares are from the offset field of 
C1J0152, filled circles from lHuang et a f| jllffi' ftl. filled squares from 
ISaracco et alJ ll997l) and stars frolii roiorgo^ki*et alJ Il995l) . 




Figure 3. Binned K band luminosity function of C1J0152 
(z=0.83). Also shown is the best fitting Schechter function. The 
solid line excludes the BCGs (one from each sub-clump)and the 
dashed line includes the BCGs. The open symbol (offset slightly 
for clarity) shows the effect of excluding the BCGs from each 
sub-clump. 



Redshift Survey jDe Propris et al.ll2003h . However, for indi- 
vidual clusters, ther e is often observed a d ip in the LF (eg. 
iBiviano etld1ll995l for the C oma cluster. lYagi et al.ll2002l 
for a sample of 10 clusters). IBiviano et all Jl995F ~suggest 
that the dip is caused by a recent episode of star-formation 
in the brightest galaxies thereby enhancing their luminosi- 
ties. However, th e presence of a dip in the H band LF of 
Coma observed bv lAndreonl (|2002) suggests that recent star- 
formation is unlikely to be the cause of the dip, since the 
H band is rather insensitive to the enha nced blue emis- 
sion accompanying star-formation. Rather, lAndreonl J2002T) 
suggests that the dip is the res ult of combining the LFs of 
several different types (see also iBingelli et al J (19881. How- 
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Figure 4. K band luminosity function of C1J1226 (z=0.89). The 
best fitting Schechter functions including (dashed line) and ex- 
cluding (solid line) the BCG are shown. 




Figure 5. K band luminosity function for C1J1415 (z=1.03). 
Schechter function fits are as for figure HI 



ever, lAndreonl i2002tl points out that such a dip is not 
seen in the near-infrar e d LF s of five clusters studied by 
iTrentham fc Mobashel lll99Sl) . nor in the H band LF of 
the Coma cluster measured by iDe Propris et "all jl99ct) . If 
the dip is the result of the combination of different LFs for 
diff erent morphological typ es, then the difference betw een 
the lDe Propris et alJ (119981) LF and the lAndreonl J2002I) LF 
suggests that the dip in the Coma LF may be dependent on 
which part of the cluster is being observed, and hence the 
particular morphological mix measured. In the present case 
no obvious dips are present in any of our LFs, and given the 
limited resolution in magnitude caused by the statistical un- 
certainties, the use of any functions more complicated than 
a Schechter function is not justified. 

Equation |H| was fit to the KLFs of all 3 cluste rs us- 
ing the maximum likelihood technique of lCashl il979T) . This 
technique fits a Schechter function to the background sub- 
tracted, binned data in the following way. The parameters 
K* and a are free parameters, and are allowed to vary be- 



tween assumed limits. For each combination of K* and a in 
the parameter space, a Schechter function is calculated, with 
the normalisation, <j>* , set such that the predicted number of 
galaxies is equal to the actual number of galaxies observed 
within the same magnitude limits. Thence the Cash statis- 
tic, C, can be computed, 



C = -2 



N 

E 



(njnei 



Inn; 



(9) 



where N is the number of bins, m is the observed num- 
ber of galaxies in the ith bin and d is the expected number 
of galaxies in the ith bin resulting from the Schechter func- 
tion. The parameters for which C is a minimum a re deemed 
to be the best fitting parameters (see lCasblll979l) . 

It was found that a could not be strongly c onstrained 
and therefore, following IDe Propris et all (^999) we fix at 
a = —0.9 which is the valu e for the KLF of the Coma clus- 
ter jDe Propris et al]ll99sl) and of the field iGardner et alJ 
1997). The characteristic magnitude, K* , was fit as a free 
parameter and 0* was determined by ensuring that the in- 
tegral of the function down to the limiting magnitude was 
equal to the total number of galaxies to the same magnitude. 

It has been noticed since early studies of cluster lu- 
minosity functions that t he brig htest cluster galaxies are 
a special class of object jPeachl IT969) and that a better 
fit to t he LF can often be made by excluding them from 
the fit iSchechterlll97fj. |Presslerlll97gl) . Therefore Schechter 
functions were fit both including the BCGs and excluding 
them. Note that in the case of C1J0152 we have excluded 
two BCGs, one from each sub-clump of the system. The 
best fitting parameters are given in tabled along with the 
X 2 probability (since the maximum likelihood statistic does 
not give a measure of the goodness of fit). The best fitting 
Schechter functions, including and excluding the BCGs, are 
also shown in figures [3] E] and 

For two of the three clusters studied here, good fits are 
found both including and excluding the BCG (see tabic |2J>. 
For C1J1226, a marginal fit is obtained excluding the BCG 
because of one high bin (at K — 20.4), but the fit is poorer 
when the BCG is included, a conclusion supported by a vi- 
sual inspection of figure 3. The values of K* from each clus- 
ter are consistent with each other, but are systematically 
fainter by ~0.2 mag when the BCGs are excluded. This dif- 
ference is of similar size to the statistical uncertainty. For 
consistency with previous measurements of Coma, in the 
following section we will refer to K* values from the fits 
excluding BCGs. 



4.2 Evolution of K* 

To evaluate the evolution in K* , models of galaxy evolution 
were ma de using the synthet i c stel lar population (SSP) li- 
braries of iBruzual fc Charlotl B2003TI . The parameters of the 
model are the star formation history, the initial mass func- 
tion and the metallicity. For all models considered here a 
single burst of star formation was assumed to occur (but 
at different ep ochs), the initial mass function was that of 
ISalpeterl 1^955;) and the metallicity was assumed to be solar 
(i.e. Z = 2%). Models were constructed for passive evolution 
in the following way. Spectral energy di s tribut ions (SEDs) 
were drawn from the lBruzual fc Charlod i2003fi libraries for 
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Table 2. Best fitting parameters of a Schechter function for each cluster with a 



A" 



Including BCGs 
(lower limit, upper limit) <f> 



Prob(x 2 ) K 



-0.9. 

Excluding BCGs 
(lower limit, upper limit) <j>* 



Prob( X 2 ) 



C10152 
C11226 
C11415 
Combined 



17.59 (17.26,17.90) 
17.79 (17.52,18.03) 
17.96 (17.63,18.26) 
17.81 (17.51,18.07) 



47.89 0.57 

57.37 0.01 

42.53 0.96 

44.69 0.07 



17.76 (17.43,18.07) 
18.04 (17.79,18.28) 
18.08 (17.75,18.38) 
17.98 (17.68,18.26) 



50.31 0.81 

61.01 0.03 

43.80 0.78 

45.59 0.16 




Figure 6. The evolution of A"*. The circles are da ta from this pa- 
per. The squares are from lDe Propris et alJ ll999f) . open symbols 
being low Lx systems and closed symbols being high Lx systems. 



different ages. A redshift of formation was chosen, corre- 
sponding to an age t = 0. Then for the adopted cosmology 
redshifts were calculated for the rest of the ages. The SEDs 
were redshifted to the appropriate redshift for each age and 
multiplied by the filter transmission curve. Thus the flux of 
the SSP could be calculated for each redshift. The fluxes 
were converted to apparent magnitudes and scaled appro- 
priately, e.g. for the evolution of K* the models were nor- 
malised to K* of the Coma cluster (see below). 

No-evolution models are made by simply choosing an 
SED of a particular age and integrating the redshifted spec- 
trum under the filter transmission curves. 

Passively evolving models with redshifts of formation 
Z{ = 1.5, Zf = 2 and Zf = 5 were made along with a no evolu- 
tion model of a lOGyr old population. The results are shown 
in figure OH The models were scaled to the Coma cluster with 
K* = 10.9 at z = 0.0231 (excluding BCGs. lDe Propris et ail 
119991) . 

It is clear that the trend seen in lDe Propris et alJ dl999T) 
is followed here and extended to higher redshifts. All values 
of K* are brighter than predicted by no evolution, and are 
generally consistent with passive evolution models. The high 
redshift points seem most consistent with zj ~ 1.5 although 
the errorbars are large enough to allow a range of Zf. A 
better fit to all the data would be obtained by normalising to 
K* at z = 0.2 (which is obtained from the mean of 3 clusters) 
rather than normalising to Coma. In that case Zf ~ 2 would 
be a better fit to the highest z data. 

Although it was not possible to constrain a for any 
individual cluster, in order to verify the choice of a = —0.9 
we have combined the counts from all three clusters into a 



general LF and fitted a Schechter function leaving a as a 
free parameter. 

The combined LF was generated by finding a weighted 
average of the counts in each bin of the individual clusters 
where the weighting factor was the reciprocal of (j>* ■ In fact 
the <j)* values for each cluster are similar, so the weighting 
has little effect. The combined LFs, excluding and including 
BCGs, with the best fitting Schechter functions are shown 
in figures and ^respectively. The best fitting values were 
K* = 18.53 and a = —0.54, when the BCGs were not in- 
cluded and q = —0.94 and K* — 17.83 when they were. 
The confidence limits, excluding BCGs, are shown in fig- 
ure |H1 where it can be seen that the constraint on a is not 
strong, ranging from a « —0.15 0.9. 

There is a large degree of degeneracy between K* and 
a as is evident from the shape of the contour in figure |H1 
Consequently when fitting a Schechter function, K * may be 
misleading if the faint end slope is poorly constrained. To 
quantify this effect the combined LF was fit with a = —0.9. 
The results are given in table|3 The result (excluding BCGs) 
is a brighter value of K* , as expected from the shape of 
figure |H1 

Consequently, if a ~ —0.54 for high redshift clusters 
then the evolution in K* may not be as strong as seen in 
figure |S| and values of 2/ ~ 2 — 5 would be more consistent 
with the K* measurements. However, if the last bin of the 
combined KLF of figure |7| is excluded then a = — 0.99 ± 
0.4. Thus, it is unclear whether the measured value of a — 
—0.54 is in artefact due to incompleteness in the faintest 
bin, and therefore we elect to use the value of a — —0.9 for 
the individual fits to clusters. 



4.3 Comparison with the Coma cluster 

We make a comparison with the Coma cluster because it 
has a similar L x (1.6xl0 45 erg s" 1 ), T x (8 keV) and mass 
(10 15 M©) as the high redshift clusters, and it is relatively 
well studied at NIR wavelengths. The K band LF of the 
Coma cluster shown in figure[7|was derived from the H band 
luminosity function of iDe Propris et alJ dl99ct) using their 
given value of H — K = 0.22 and to a magnitude limit for 
which all galaxies had spectroscopic membership confirma- 
tion. Note that the two LFs in figurc[7]have not been scaled 
in any way; the high redshift clusters are of similar richness 
as Coma (although the comparison is not exact because the 
high redshift LFs were obtained from a larger fraction of 
t he cluster area than for Coma; see below). The Coma LF 
ofjAjifL]e£n_^£l^lk2|j [2000l) covers a smaller area than that of 
IDe Propris et al.l II1998T) but is in good agreement (except for 
a dip in one bin). For the high redshift clusters, the conver- 
sion from apparent to absolute magnitude was made using 
a fc-correction calculated for a lOGyr old synthetic stellar 
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Figure 7. The combined luminosity function of all three clusters 
((z)=0.9) shown by solid circles, excluding BCGs. Its best fitting 
Schechter function is the solid line, whilst the dotted line shows 
the expected evolution at the redshift of Coma assuming passive 
evolution and Zf = 2 . The number of galaxies in each bin is a 
weig hted average. The Coma K band luminosity function derived 
from lDe Propris et alJ Jl99Sfl is shown by open squares. Note that 
the two LFs are plotted directly as observed, with no scaling. See 
text for details 



population from the libraries of iBruzual fe Char rod (2003). 
The typical size of the fc-correction was « —0.7 mag. 

The offset in absolute magnitude of K sal. 2 mag be- 
tween the Coma KLF and the combined high redshift KLF 
indicates the degree of evolution. The same degree of evolu- 
tion is also seen in figure |S| as the difference between the ob- 
served values of K* and the no-evolution predictions. Evolv- 
ing the fit to the high redshift KLF, assuming a passively 
evolving model with Zf = 2, gives a surprisingly good de- 
scription of the observed Coma KLF (the dotted lines in 
figures |7| & |UJ . Thus differences between the high redshift 
and low redshift systems, as shown by the bright end of 
the LF studied here, can be accounted for by purely passive 
evolution. There is no obvious change in the shape of the 
LFs. 



Figure 8. Confidence limits at 68 % on the best fitting param- 
eters for the combined luminosity function when BCGs are ex- 
cluded. 




Figure 9. The combined luminosity function of three clusters 
including BCGs. Symbols are as in figure 171 



4.4 Integrated light functions. 

An interesting and complementary measure of galaxy evo- 
lution within clusters is provided by their integrated light 
function. This is calculated from the LF using 

fe 

j(<L fe ) = ^± k = l,...,n (10) 

i=l 

where iV; is the number of galaxies in the i th bin of the 
LF and L; is the luminosity of the bin, and n is the total 
number of bins. For the BCGs and bins in which there is 
only one galaxy the individual contribution was calculated 
using the total luminosity of the galaxy and not the bin 
centre. 



Integrated light functions were calculated for all three 
clusters and for the Coma cluster. These are shown in fig- 
ures ^| and ^3 f° r data excluding and including BCGs. In 
order to compare each cluster the luminosity has been nor- 
malised to I*. Thus the effects of any pure luminosity evo- 
lution (such as passive evolution), which would alter the lu- 
minosity by a constant factor, are eliminated from the plot 
and any variation between the shapes is due to some other 
evolutionary effect. Note also that in order to make this com- 
parison the integrations of the light functions must reach the 
same depth below L* . As this was not the case the counts 
for the faintest bin for C1J1415 and the two faintest bins 
for C1J0152 were calculated from the best fitting Schechter 
functions and thus these are not shown in the plots. The 
contribution from these bins is only 5-10% of the light inte- 
grated to our limit of L/L* — 0.05. Integrating the Schechter 
function to luminosities fainter than this limit would add 
only a further « 1%. 

A comparison of the shapes of the integrated light func- 
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Figure 10. Integrated light functions for all three clusters com- 
pared to Coma. Brightest cluster galaxies have been excluded. 



reach radii of 30%-40% of the virial radii (as measured from 
the X-ray properties and provided by Ben Maughan, pri- 
vate communication) , whereas the observation of the Coma 
cluster covers only 0.16 of the virial radius (r2oo = 2 . 3 Mp c 
for Hp = 70km s " 1 Mpc" 1 . iReiprich fc Bohringerl [2002). 
iDriver et alJ l|l99Sfl show that the dwarf to giant galaxy ra- 
tio within rich clusters is strongly anti-correlated with the 
mean projected density outside the cores of clusters. Our ob- 
servations donot, however, reach dwarf galaxy luminosities. 
lAndreor] |2001) confirms that the fraction of dwarf galaxies 
increases in the outer parts of a z — 0.3 cluster, AC118, 
and furthermore shows that the faint end slope of the KLF 
recovered from the outer parts of the cluster is considerably 
steeper than the faint end slope of the KLF recovered from 
the core. Consequently comparisons of observations cover- 
ing different fractions of the virial radius must be made with 
these effects in mind. However, our concern is with the bright 
end of the LF in the inner regions of the clusters, where these 
effects are unlikely to produce large differences. 




Figure 11. Integrated light functions as in figure IT0l but with 
BCGs included. 

tions for the four clusters excluding their BCGs, shown in 
figure ITUl shows that they are all consistent. If the BCGs 
are included in the integrated light function, as in figure ITT1 
the light functions are generally consistent within errors but 
the scatter is greater. 

If merging plays an important role in the evolution of 
massive galaxies since z « 0.9 then it would be expected 
that there would be a higher fraction of total cluster light 
in bright galaxies at low redshift compared to high redshift. 
Little evidence is seen for this in either figure 1101 or 1111 If 
BCGs are included then C1J1226 is seen to have a slightly 
higher fraction of light in its bright galaxies than Coma and 
C1J0152 has a lower fraction. The lack of any systematic 
difference between the integrated light functions of the high 
redshift clusters and Coma is suggestive that merging does 
not play a dominant role in the evolution of massive early- 
type galaxies since z w 0.9. 

One caveat is that interpretation of the integrated light 
functions is complicated by the different areas observed for 
each cluster. The observations of the high redshift clusters 



5 THE BRIGHTEST CLUSTER GALAXIES. 

The uneasy fitting of some BCGs on the Schechter function 
may be symptomatic of a different evolutionary process to 
the majority of galaxies within a cluster. We now examine 
this evolution. 

The K band Hubble diagram is shown in figure 1121 
The models are the same as described above, only now 
they are normalised to the low redshift [z < 0.1), high Lx 
(> 1.9 x 10 44 ergs s " 1 0.3 - 3.5 keV) BCG photometry of 
iBrough et alJ i2002T) . The high Lx clusters from figure 1 
of iBrough et alJ J2002D are shown as the open circles. The 
BCGs in the current study, shown as solid points, extend 
to higher redshifts. Note that there are two BCGs shown 
for C1J0152, one for each sub-clump. The magnitudes of 
the BCGs were measured in a metri c aperture of radius , 
r m = 12.5ft _1 kpc for comparison with IBrough et ail |2002). 
The BCG in C1J1226 overlaps with at least one fainter 
galaxy observed in projection; this galaxy has been sepa- 
rated by Sextractor. We note that the total magnitude of the 
C1J1226 BCG measured here (7^=15.34 mag) is consistent 
with the isophotal value of if=15.5 measured by Cagnoni et 
al (2001). 

It can be seen from figure lT^l that the BCGs in C1J1226 
and C1J1415 are somewhat brighter than no-evolution pre- 
dictions but in good agreement with passive evolution mod- 
els. The BCGs of C1J0152 are fainter relative to the models 
than any other BCGs in figure [T2l at z > 0.4, perhaps indica- 
tive that the BCGs are not ful ly formed. Indeed as C 1J0152 
is probably a merging system jMaughan et a l. 2003) it may 
be supposed that if the BCGs merge, and no further star- 
formation is triggered, they will produce a brighter BCG 
more consistent with the other high redshift BCGs. 

In contrast C1J1226 has a relaxed X -ray morphology 
jMaughan et all2005trEbeling et ai]l200ll) . and is one of the 
brightest BCGs at z > 0.4 relative to the models. These 
results are suggestive that the evolution of BCGs at high 
redshifts may be related to the dynamical history of the 
cluster they inhabit, although confirmation awaits a larger 
sample. 

Table lists the contributions from the BCG to the to- 
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Figure 12. A" band Hubble diagram for BCGs, So l id poi nts are 
from this paper, open points are from lBrough et al.l J2002I) . Mag- 
nitudes are measured within a 12.5h~ 1 kpc radius aperture. 



Coma with no need to hypothesise any significant mergers 
in their future evolution. 

Why are the BCGs in C1J0152 significantly fainter than 
the BCG in C1J1226? We will test the assumption that the 
stellar mass of the BCG is related to the mass of its host clus- 
ter. We assume the ratio of masses of the BCGs in C1J0152 
and C1J1226 is proportional to the ratio of the total cluster 
masses times the ratio of the fractions of galaxy light in the 
BCGs, assuming the mass fraction in galaxies is the same 
in each cluster. Then the combined effects of the lower to- 
tal mass of each subclump of C1J0152 (by a factor of ~0.39 
compared to C1J1226) and the lower fraction of cluster K- 
band light in the BCGs (by a factor of «0.69) are consistent 
with the difference in the luminosities of the BCGs (by 1.57 
mag or a factor of 0.24). This suggests that the BCG stellar 
masses at high red shift may be lower in l ower mass host clus- 
ters , as inferred bylBrough et al.l d2002l) . lBurke et al.1 (l200(t) 
and lNelson et alJ (1200211 . 

We note that if the BCGs in C1J0152 were to merge, 
as was hypothesised to explain their positions in the Hub- 
ble diagram, the fraction of light in the BCG would remain 
roughly similar. 



Table 3. Fraction of the cluster i^-band light in the BCGs. The 
values for the BCGs in C1J0152 are the fractions in each sub- 
clump independently. 

Cluster Fraction Radius/ r v j r ^"bcg 



C1J0152A 0.10 0.29 16.87 

C1J0152B 0.12 0.34 16.62 

C1J1226 0.16 0.41 15.34 

C1J1415 0.11 0.43 16.27 

Coma 0.06 0.37* 8.86 



* extrapolated from r=0.17r v ; r , see text. 



tal K band light for each cluster. The total luminosity in 
each cluster was calculated from the LFs, as for the inte- 
grated light functions, to a depth of L/L* w 0.06. The sub- 
clumps in C1J0152 have been treated independently, since 
each has a BCG. In the case of Coma a correction was made 
to account for the smaller fraction of the virial radius ob- 
served ;_riiis_was_cateubated using the number density pro- 
file of iKent fc Gunnl (1982) to determine the ratio of the 
number of galaxies within the observed area to that within 
0.37r v ir, the average area observed for the high redshift sam- 
ple. The correction was a factor of 1.96. The luminosity of 
the Coma BCG, NGC4874, was calculated using the 2MASS 
total magnitude. Note that there are two galaxies of simi- 
lar brightness in Coma, NGC4874 and NGC4884. We have 
taken NGC4874 as the BCG here as it resides at the centre 
of the X-ray emission. 

The results show that the fraction of light in the BCGs 
at high redshift is greater than, or equal to, the fraction in 
Coma. This is still true if both brightest galaxies in Coma 
are included. 

There are uncertainties in the corrections applied in or- 
der to compare the clusters, which were originally measured 
to slightly different depths and covered different areas, al- 
though the variations in radius would produce only small 
variations in total K band light (~20% using the Coma light 
profile). A physical interpretation is that the high redshift 
BCGs are already of similar or greater mass to the BCG in 



6 DISCUSSION AND CONCLUSIONS 

The evolution of the galaxy populations of three high red- 
shift clusters of galaxies has been studied. The bulk evolu- 
tion of the galaxies, as characterised by K* , is found to be 
consistent with passive evolution with a redshift of forma- 
tion Zf ~ 1.5-2. Further evidence for passive evolution is 
seen in the similarity of the shape of the high-redshift LF 
with that of Coma, and in the consistent shapes of th e in- 
tegrated light functions. iTrentham fc Mobasherl dl998h also 
reached similar conclusions about the evolution of the shape 
of the isT-band LF, albeit with poorer accuracy. 

Purely passive evolution of early-type galaxies is con- 
sistent with several other studies i ncluding the evolution o f 
the K band luminosity function dDe Propris et aDll999T) . 
evoluti on of the fundamental p lane in terms of mass-light 
ratios dvan Dokkum et alJll99cth and stud ies of the scatter 
of the colour-magnit ude relation (see e.g. lEllis et aj]ll997L 
IStanford etaflll99Sl) . 

Our conclusions are different to those of iBarger et alJ 

( 1998), who found no evidence of evolution in for clus- 
ters at 0.31<z<0.56 assuming a qo=0.5 cosmology. However, 
the cosmological dependence is such that at these redshifts 
the qo=0.5 NE prediction is very similar to the (f2 m =0.3, 
Q A =0.7) Zf = 2 passive evo lution prediction (see e.g. figure 
8 of iDe Propris et aljll999T) . Thus the results of Barger et 
al. are in agreement with those found here. 

When discussing formation it is important to distin- 
guish between the epoch at which the stars in the galax- 
ies were formed and the epoch at which the galaxies were 
assembled. The studies of the fundamental plane and the 
colour-magnitude relation refer to the epoch of star forma- 
tion. If merging were a dissipationless process then it would 
be possible to have no extra star formation as a result of a 
merger and thus the age of the stars within a galaxy could 
be older than the age of galaxy assembly. A study of the 
cluster of galaxies MS 10 54-0 3 at z = 0.83 is pre s ented by 
Ivan Dokkum et alJ dl999l) and Ivan Dokkum etafl (1200(f) m 
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which there is observed a high fraction of merging red galax- 
ies. Very little star formation is seen in the merging galaxies 
constituting evidence that the galaxies are in fact somewhat 
younger than the stars that reside within them. 

Is such merging reflected in the evolution of the LF? 
The K band luminosity of a galaxy is very nearly indepen- 
dent of star-formation, but reflects the mass of the old stars 
within the galaxy. Thus K mag nitudes are a good me asure 
of the stellar mass of a galaxy. iDiaferio et al.l i200ll) give 
predictions of the evolution of the KLF from semi-analytic 
models of dissipationless, hierarchical structure formation. 
The models show that there is very little evolution of the 
number of massive galaxies in clusters since z = 0.8. The 
galaxies are assembled at high redshift and evolve passively 
with little subsequent merging after z = 0.8. A detailed 
comparison with these models is however not possible be- 
cause the luminosity evolution predictions in the models are 
not s ufficiently accurate for massive galaxies fDiaferio et alJ 
120011) . 

It is clear from figure [5] that the bulk of galaxies in our 
sample were brighter in the past than predicted from no- 
evolution models, by AK ~ -1.2 mag at z = 0.9. A direct 
comparison of the high-redshift K band LF with that of 
Coma suggests that the two are very similar in shape and 
the differences in may be reconciled by pure luminos- 
ity evolution. The fading with time by AK « -1.2 mag is 
consistent with passive evolution from a formation epoch 
Zf ~ 2. In the monolithic collap se picture this would be ex- 
pected naturally. The models of lDiaferio et al.l i200l|) show 
that this is also predicted for massive cluster galaxies in a 
hierarchical scenario since most merging takes place early 
on the history of the cluster. 

In a merging model passive evolutionary processes will 
still be present, and thus K* would still appear brighter than 
no-evolution predictions. A probe of 'extrapassive' processes 
is the shape of the LF. It is found that a is consistent with 
that of Coma although it is poorly constrained here. Perhaps 
a stronger test for extrapassive processes is the shape of the 
integrated light function. The lack of any major changes 
seen in figure tTUl is suggestive that passive evolution alone is 
responsible for the evolution measured in K* . We conclude 
that the luminosity evolution of bright galaxies in massive 
clusters is consistent with pure passive evolution, but note 
that this may be consistent with hierarchical models if most 
merging takes place at high redshifts. 

The mild if-band evolution of luminous field galax- 
ies to z=l found in recent surveys is also consistent with 
luminosity evolution, alth ough the amount of evolution 
AM* K = -0.54± 0.12 mag JPozzetti et alJl2003ft or AM k = 
-0. 7 ± 0.3 mag jFeulner et al.ll2003t iDrorv et aT]|200l see 
also llm et al][2 002'l is le ss than that observed her e (AM^ = 
-1.2 ± 0.3 mag) or bv iDe Propris et alJ lll999D (AMJ = 
—0.90 ± 0.25) in massive clusters. The recent results on the 
field evolution indicate little change i n <j>* to z=l, and are in 
contra st to the conclusions d r awn b y lKauffmann &: Chariot! 
( 1998) and Kauffm ann et all lll996f) who found evidence for 
a deficit of massive galaxies in the field at z ~ 1. 

The difference in luminosity evolution between the field 
and massive clusters (both at z«0.9) is in the opposite sense 
to the environmental dependence of the formation epoch 
in hierarchical models <|Kauffmann et alJll993t iBaugh et alJ 
ll996l.lDiaferio et alJl200ll) . In these models, the assembly of 



an early-type galaxy, and the formation of its stars, occurs 
at an earlier epoch in a region of high overdensity than in 
the field. The degree of luminosity evolution between z=l 
and z=0 is predicted, therefore, to be higher in early-type 
field galaxies than in massive clusters, since z=l is closer to 
their epoch of formation. If the high redshift measurements 
are not in error (they can be checked via spectroscopic sam- 
ples), then more detailed studies will be required to separate 
the if-band evolution of galaxies of different spectral types, 
since the cluster and field samples contain a different mix of 
galaxies and thus a different mix of star formation histories. 

Note that although there is a strong case that the ob- 
served evolution in K* may be attributed almost entirely to 
passive evolutionary processes, the interpretation of this re- 
sult as being due to a redshift of formation Zf = 1.5 — 2 
is less secure. The models used are for a single burst of 
star-formation for a stellar population with a Salpeter ini- 
tial mass function having a solar metallicity. Age-redshift 
relations are calculated for an assumed flat cosmology with 
H = 70km s" 1 Mpc" 1 and Q M = 0.3 and fi A = 0.7. All 
of these assumptions affect the resultant redshift of forma- 
tion and so there is clearly some slack in the interpretation 
of the evolution. For example if the models if figure |S| are 
normalised to the high Lx point at z — 0.2, rather than 
Coma, the data are more consistent with Zf = 2-5. Future 
work based on the colours and morphologies of the member 
galaxies will provide stronger constraints on the epoch of 
formation. 

The evolution o f BCGs is consistent with that found by 
iBrough et alJ {2002) for high Lx clusters. Figure H2lexhibits 
a degree of scatter in the evolution of BCGs at high redshift 
with some BCGs being consistent with no-evolution pre- 
dictions, of which C1J0152 is an example, and others being 
consistent with passive evolution. This can be interpreted as 
showing that some BCGs are fully formed at high redshift 
e.g. C1J1226, whereas others would need to undergo merging 
between z — 1 and the present to reconcile them with local 
BCGs. 

This is supported by the ratios of K band light in the 
BCGs to total cluster light. The ratio in Coma is not larger 
than the ratio in the high redshift clusters. Thus excess 
brightening of the BCGs with time, relative to the other 
cluster galaxies, as would be expected from such processes 
as cannibalism, is not observed, suggesting that the BCGs, 
at least of C1J1226 and C1J1415, are already fully formed. 
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